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1. Introduction 



Trapped ions have featured prominently in the study of elementary quantum systems over a period of 
many years. Some of the earliest experiments demonstrated the phenomenon of electron shelving or 
quantum jumps [1-3]. Subsequently, with improved methods of laser cooling [4-6], and through resolved- 
sideband Raman cooling [7,8] in particular, it became possible to prepare an ion in the ground state of a 
harmonic trap. Thus, a broad range of elementary quantum physics was opened up to experimental study 
^ , [9], including the spectroscopy [10-12] and photon statistics [13-15] of single-atom resonance fluorescence; 

' quantized Rabi oscillation and engineered quantum states of motion [16-18] (through a "phonon" analogy 

^ , of the Jaynes-Cummings model); and quantum decoherence, specifically, of superpositions of motional 

states [19,20]. Today, trapped ion systems continue to be developed as potential building blocks of 
quantum information processors [21-23]. 

Theoretical work on the quantized motion of trapped ions has mostly considered small displacements 
from equilibrium — the so-called Lamb-Dicke regime [24]. It is common for the ion to be manipulated in 
a standing- wave optical potential [25,26], in which case the potential is expanded to lowest order in the 
, Lamb-Dicke parameter 

27r 



> 
o 



\ where rj is the ratio of the ion ground-state position uncertainty Ax — ^jTil2mujT and the wavelength 

' of the standing wave A (multiplied by 27r); m is the ion mass and ujt the frequency. With the restriction 

to the Lamb-Dicke regime, in many situations only one, two, or possibly a few quanta of excitation are 
^ ■ considered [16,17,27,28]. 

rS \ The assumption -q <^ 1 is a feature of the semiclassical treatment of atomic motion in laser light as 

■ well, where it underlies the diffusion or Fokker-Planck equation models [4,29-31]. These models assume 

that significant change to the particle (ion) momentum is built up from very many individual momentum 
"kicks" (each of size fik = h/X) with the wavepacket well-localized on the scale of A; thus, in a harmonic 
trap it is assumed that 

Ap/hk = {2kAx)-'^ = {2t])-^ > 1. (2) 

Such treatments also adiabatically eliminate the internal degrees of freedom, requiring the distance 
traveled in an atomic lifetime to be much less than A. This imposes the additional restriction 

iujT/l)xo < A, (3) 

where xq is the amplitude of oscillation in the trap. Since, reasonably, xq > Ax — riX/2TT, from Eqs. (2) 
and (3), the recoil energy is constrained by 

= ■q^uoT l < 27r7?. (4) 

Ti'y 



In this paper wc; study the quantum motion of a trapped ion in laser light under eonditions where none 
of the usual limitations and restrictions hold. Our aim is to exhibit the trapped ion system as an example 
of a quantum stochastic process in its fully elaborated form. Specifically, we consider an ion driven on a 
two-state electronic resonance by a standing- wave laser field — i.e., the heating of an ion through resonance 
fluorescence in a standing wave. Motion in one dimension is considered under the assumption of strong 
transverse confinement. The following features are notable in relation to the limitations and restrictions 
above: 

(i) heating rather than cooling of the ion takes place; thus, beginning in the ground state of the trap, 
the amplitude of oscillation grows to eventually extend over many periods of the standing wave, 
where energies (depending on r]) as high as IO^Ucot — 10^ fiUT are reached; thus, quantum motion 
well outside the Lamb-Dicke regime is considered; 

(ii) we focus upon the case 

wt/7 = 1, (5) 

so inequality (3) is violated at relatively small amplitudes of oscillation; a crucial interplay develops 
between the coherent evolution of the center-of-mass and internal degrees of freedom; even when a 
diffusion paradigm fits (for sufficiently small rj), the motion is not understandable as diffusion in a 

prescribed harmonic plus optical potential; 

(iii) strong laser excitation, with Rabi frequency 

= 2, (6) 

is considered, but due to the interplay just noted the dressed-atom approach to atomic motion in 
laser light [32] cannot be used; modulation of the Rabi frequency by the harmonic center-of-mass 
motion plays a central role; 

(iv) Lamb-Dicke parameters ranging from 77 = 0.2 to = 3.0 are considered, thus demonstrating the 

transition from a quasi-classical dynamic, for 77 <C 1, to a manifestly quantum, stochastic dynamic 
when ?7 > 1; this transition is similar to that from weak- to strong-coupling in the treatment of 
quantum noise in cavity QED [33]. 

The reported investigation makes use of full quantum trajectory simulations [34], quantizing both 

the internal and the center-of-mass degrees of freedom. As an aid to understanding, approximate, 
semi-quantum trajectory simulations are explored as well, where the center-of-mass motion is treated 
classically apart from the inclusion of momentum "kicks" coordinated with the fluorescence. There have 
been previous quantum trajectory (Monte Carlo wave- function [35]) simulations of atomic motion in laser 
light [36-38], but all, to our knowledge, were carried out for conditions satisfying the restrictions of 
Eqs. (2)-(4). 

We begin in Sec. 2 with an overview, presenting a summary of results in the form of heating curves 
(mean energy of the ion as a function of time) calculated for a series of Lamb-Dicke parameters ranging 
between 77 = 0.2 and 77 — 3.0. Results calculated from semi-quantum and quantum trajectories are 
compared and a number of differences identified for explanation in subsequent sections. TIk^ equations 
underlying the Monte-Carlo simulations are presented in Sec. 3, where we demonstrate a surprising 
feature of the individual realizations of the stochastic process; a stepwise evolution of the amplitude of 
oscillation of the ion is observed, pointing to the (existence of a series of mctastable amplitudes on which 
the ion heating almost stops. In Sec. 4 the origin of this evolution is identified as frequency modulation 
of the internal state Rabi oscillation. The metastable amplitudes of oscillation are identified with zeros of 
the Jq Besscl hmction and characterized by the mean waiting time for fluorescent scattering plotted as a 
function of the amplitude of oscillation in the trap. A phase-space diffusion model is derived that recovers 
the heating curves of Sec. 2 in the limit of small Lamb-Dicke parameters, and quantitative differences 
between the semi-quantum and quantum diffusion are explained. Finally, in See. 5, Lamb-Dicke 
parameters of order unity are considered, where differences between a diffusion and a jump process 
are revealed. Here semi-quantum and quantum trajectory models are qualitatively different in their 
predictions. An observed quantum suppression of the ion heating rate is explained by the delocalization 
of the center-of-mass wavepacket across the standing wave. 



2. Quantum and semi-quantum heating rates 

We consider an ion of mass m trapped iiarmonically in one dimension with trap frequency lot- A classical 
laser field resonantly excites a closed two-state electronic transition of the ion, with decay rate 7 and Rabi 
frequency f2; parameters chosen as in Eqs. (5) and (6). The laser field forms a standing wave along the 
axis of the ion motion (wavelength A) and the equilibrium of the potential is located at an anti-node of the 
standing wave. 

Figure 1 shows the growth of the mean-squared amplitude of oscillation of the ion — measured in wave- 
lengths [see Eqs. (20) and (21)] — as a function of time for various Lamb-Dicke parameters r]. Two separate 
models were used to compute these results as ensemble averages over Monte-Carlo simulations of the ion 
heating through fluorescence. For the results of Fig. 1(a), the ion center-of-mass position and momentum 
are treated as classical variables, but with a momentum "kick" added when a laser photon is scattered; the 
scattering is simulated as a quantum trajectory for a two-state system with a definite, though oscillating, 
position within the standing wave [34-37,39]. For the results of Fig. 1(b), full quantum trajectory 
simulations were carried out with both the center-of-mass and internal degrees of freedom quantized. The 
following points are of note and will be elaborated upon in the following sections: 

(i) for small Lamb-Dicke parameters both sets of results approach a limit where the heating curve 
becomes independent of r/, apart from a factor rf that may be absorbed in the scaling of time; in 
this limit the heating is well-described by a diffusion in phase-space; 

(ii) semi-quantum and quantum trajectory models disagree quantitatively in the diffusion limit 
[compare the thinest curves in Figs. 1(a) and (b)]; 

(iii) larger Lamb-Dicke parameters yield an 77-dependence over and above the scaling of time by 77^; thus, 
differences between a quantum jump and diffusion process become explicitly apparent; 

(iv) for large Lamb-Dicke parameters the two sets of curves are qualitatively different; the semi-quantum 
trajectories [frame (a)] show a heating rate that increases monotonically with increasing -q, while the 
full quantum trajectory results [frame (b)] eventually show a dramatic reduction of the heating rate. 

Before taking up these four themes, we first present the mathematical details of our two models in 
Sec. 3. Then we explore individual realizations of the quantum stochastic heating of the ion. These 
exhibit a very interesting evolution of the amplitude of oscillation of the ion in the trap, an underlying 
structure to the stochastic dynamics quite unanticipated on the basis of the ensemble averages displayed in 
Fig. 1. 




Fig. 1. Heating of a harmonically trapped ion through resonance fluorescence in a standing- 
wave field. Plots of the mean squared amplitude of oscillation as a function of time from: (a) 
semi-quantum trajectory simulations — Lamb-Dicke parameters rj = 0.2, 0.4, 0.8, 1.4, 2.2, and 
2.6 (thinest line to thickest line): (b) quantum trajectory simulations — Lamb-Dicke parameters 
)7 = 0.8, 1.4, 2.2, 2.6, 2.8, and 3.0 (thinest line to thickest line). 



3. Monte-Carlo simulations 



3.1 Semi-quantum trajectories: 



Semi-quantum trajectories treat the center-of-mass motion of the ion classicaUy. Only the electronic states 
of the ion are quantized. We represent the center-of-mass state by the point in phase-space 



a{t) = a{t)e 



muiT 
2h 



x{t) + i 



Pit) 



mu)T 



(7) 



where x{t) denotes the oscillator position and p{t) its momentum. The trajectory equations are 
formulated in a rotating frame, with a{t) denoting the complex amplitude of the oscillation. Following 
the standard quantum trajectory treatment of resonance fluorescence (see Sec. 8.2 of Ref. [34]), the 
conditional internal state of the ion is expanded in the interaction picture as 



|V'REc(t))=c(-)(t)|-)+cW(t)|+) 



;(+)/ 



(8) 



where |— ) and |-|-) are, respectively, the ground and excited electronic states, and the label REC refers 
to the conditioning of the state on the scattering record prior to time t; ket vectors and Schrodinger 
amplitudes with an overbar, as in Eq. (8) , are unnormalized. Under the assumed resonant excitation, 
in between photon scattering events, the (unnormalized) Schrodinger amplitudes obey the equations of 
motion 



dt 
dcM 

dt 



= — {p,/2) cos\ri{ae 



+ a'e 



)] 



= (0/2) cos[r?(ae-''"-^* + a*e''"-^*)] c' 



(9) 



where, importantly, the Rabi frequency f2 is modulated by the harmonic motion of the ion through the 
standing wave — using Eqs. (1) and (7), cos[kx{t)] = cos{ri[a{t) + a*{t)]}. Photon scattering is accounted 
for through quantum jumps, which interrupt this coherent evolution. The jumps occur at random times, 
with instantaneous rate 

|c-(+)(t)p 



Each jump returns the ion to the ground electronic state. 



(10) 



(11a) 



and at the same time we add a momentum "kick" to the center-of-mass. The magnitude of this "kick" is 

determined by projecting the recoil for photon scattering in a randomly chosen direction (9, (jy) onto the 
line of the ion motion; thus, for photon scattering at time t, we implement a momentum "kick" in the 
rotating frame through the semi-quantum jump 



a + irj sin 6 cos (l>e 



iujTt 



(lib) 



which follows from the magnitude of the recoil fik and Eqs. (1) and (7). The angle (p is uniformly 
distributed between and 27r, while is distributed according to the dipole radiation pattern (oriented 
perpendicular to the ion motion). 

A Monte-Carlo simulation must decide at each time step whether to propagate the internal state 
according to Eq. (9), leaving a unchanged, or to implement the quantum and semi-quantum jumps, 
Eqs. (11a) and (lib); the branching ratio is given by the jump probability i^jumpAt. When a jump 
(photon scattering) does occur, two random numbers determine the angles 9 and (p. The heating of the 
ion results from the accumulated momentum "kicks" (lib). The mechanism is therefore, in principle, 
straightforward. It is important to note, however, that the "kick" rate is not a fixed fmiction of either 
space or time, but itself evolves stochastically according to Eqs. (9) and (10). This gives rise to an 
unanticipated complexity in the ion motion, as illustrated in Sec. 3.4. 



3.2 Quantum trajectories: 



The full quantum trajectory treatment of the ion heating is formulated as a natural extension of Eqs. (8)- 
(11b). The principle difference is that Schrodinger amplitudes c^~\t) and c^^\t) are replaced by center- 
of-mass ket vectors, with the expansion of the total quantum state, quantized internal and center-of-mass 
degrees of freedom, 

IV^REcW) = IVS^E^cC*))!-) + \i^^n^cm+)- (12) 

The coupled equations of motion (9) for Schrodinger amplitudes are replaced by coupled equations of 
motion for the (unnormalized) center-of-mass kets, 

= -(0/2) cos[^(ae-^-* + fite^-*)] 1^;^^), 

= (0/2)cos[r?(ae— * + ate^-*)] \^^-^^) - (13) 

where d and are oscillator annihilation and creation operators, and replace the classical phase-space 
variables a and a*. In line with the rotating frame of reference used in Eqs. (7) and (9), the interaction 
picture of the harmonic oscillator has been adopted in Eqs. (12) and (13). The jump (photon scattering) 
rate is calculated from the center-of-mass ket norms in a straightforward generalization of Eq. (10), 

R (t) = (^^cmnJcit)) (14) 

and (11a) and (lib) are combined in the quantum jump 

IV^REc) ^ ^(irysin^cos.^e''^-*) |V^^t;c)|-)' (15) 

where -D(^) = exp(^a''' — ^*a) is the displacement operator. 

Monte-Carlo simulations may be carried out on the basis of Eqs. (13)-(15) in essentially the same 
manner as before. The numerical demands are significantly higher, though, due to the much larger 
number basis states involved. The two center-of-mass kets are expanded in the energy representation of 
the harmonic oscillator, and we have, in fact, carried out simulations where the energy of the center-of- 
mass motion rises to be of the order of W^Tilot — IO^Tiuit- This is achieved by carrying out the 
calculations in a local frame that tracks the average complex amplitude of the oscillation, (a(t))R,EC, 
reached at the conclusion of a quantum jump. This strategy also circumvents the need to explicitly 
implement the displacement in Eq. (15). 

3.3 Quantum trajectories in the local frame: 

We absorb the mean displacement of the center-of-mass kets from one quantum jump to another into a 
change of the frame of reference. Thus, we work with the displaced center-of-mass kets, 

lCii(*)) = D[-Htm^K^cit)), (16) 

where D[—a{t)] displaces the center-of-mass state to a "local frame". The local frame changes in time, 

discontinuously, at the times of the quantum jumps. When it has been decided that a photon scattering 
is to take place and the associated quantum jump (15) is to be executed, a displacement of the post-jump 
kets is made to first cancel the displacement of the momentum "kick" [Eq. (15)], and second to move the 
mean complex amplitude of the oscillator back to zero — i.e., to off-set any change in the mean brought 
about by the Schrodinger evolution during the interval that has elapsed since the last quantum jump. 
Then in place of Eqs. (13) we have equations of motion for the ket vectors in the local frame, 

= -(0/2) cos{v [{a + a)e— + {a* + at)e--*] } l^if+i,), 

= (0/2) cos{,7 [(a + a)e-^--* + {a* + at)e--*] } li^Uj - ^li^tii), 

(17) 



with jump rate 



(^•LiliOk'Lal(O) 



T,(+) 



(18) 



The quantum jump is now executed by setting the mean oscillator amplitude with respect to the local 
frame to zero (and setting the ion in the ground electronic state), with 



^(-(a)local)|ViLll), 



0, 



and moving the location of the local relative to the global frame, with 

a ^ a + (a)iocai + iv sin 6 cos 0e"^'^*. 



(19a) 



(19b) 



In these expression (a)iocai is the conditional oscillator amplitude expectation calculated in the pre-jump 
local frame. 

3.4 Sample results: 

Figure 2 presents examples of the simulated trajectories for four different values of the Lamb-Dicke 
parameter, increasing from r] = 0.2 in frame (a) to 77 = 2.2 in frame (d). We plot the amplitude of 
oscillation, measured in optical wavelengths, against time in atomic lifetimes scaled by 7]'^ (see Sec. 4.2). 
Explicitly, for semi-quantum trajectories, using Eqs. (1) and (7), the amplitude is given by 



A V mujT 



a\ = ^\a\ 



while for quantum trajectories 

A = 



\J (a'''a)REC = '^\J ((<5* + a)){a + a))iocai- 



(20) 



(21) 



There are two principal points of interest to be noted in the figure. First, frames (a) and (b) shown a clear 
step-like evolution, which is at least qualitatively similar in the semi-quantum and quantum simulations. 
Long periods of almost no heating at all (mctastablc amplitudes) arc present, interspersed with periods 
of much more rapid heating, thus making up the steps. Quantitatively, periods of essentially no heating 
are observed to last as long as 10^ lifetimes, while the metastable amplitudes are spaced approximately 
A/2 apart — at A ~ 0.38,0.88, 1.38, etc.. Second, the similarity between semi-quantum and quantum 
trajectories seems to be breaking down in fames (c) and (d). The quantum trajectory shows no clear steps 
in frame (d), while the stepping of the semi-quantum trajectory no longer keeps to the metastable 
amplitudes identified before. The task of the following sections is to explain these features along with 
points (i)-(iv) of Sec. 2. 
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Fig. 2. Sample quantum (thick curves) and semi-quantum (tliin ciuvcs) trajectory simulations of the amplitude of 
oscillation in the trap; for Lamb-Dicke parameters rj = 0.2 (a), 0.8 (b), 1.4 (c), and 2.2 (d). 



4. Mean waiting times cind the diffusion limit 



4-1 Semi- quantum waiting times: 

An explanation of the step-like evolution of Figs. 2(a) and (b) is provided by Fig. 3. Here, in the upper 
series of pictures, we plot the induced Rabi oscillation and simulated photon scattering sequence for four 
different amplitudes of oscillation in the trap. For simplicity, only Semi-quantum trajectories are consid- 
ered in this section. Frames (a) and (c) show the ion being fully excited, while (b) and (d) show much 
lower levels of excitation. The corresponding numbers of scattered photons are high — frames (a) and (c)- 
and much lower — frames (b) and (d). Evidently, the response of the ion to the modulated driving field 
[Eqs. (9)], 



(0/2) cos[kx{t)] = (0/2) cos{27r^cos[a;Tt - arg(a)]}, 



(22) 



depends strongly on the amplitude of oscillation in the trap. The observation is readily understood from 
the time-averaged Rabi frequency, OJo(27r^), which is zero at the zeros of the Jq Bessel function. The 
first two of these zeros occur at the amplitudes A = 0.38 and A = 0.88 of frames (b) and (d) in the figure, 
while amplitudes A — 0.0 and A = 0.6 [frames (a) and (c)] correspond to local maxima of the Bessel 
function. Thus, the metastable amplitudes are located by the zeros of Jo{2ttA), where the ion is only 
weakly excited and the fluorescence rate drops. Of course, some photon scattering occurs at these zeros, 
since the actual response is to the modulated driving and not to its time average; the amount decreases as 
A increases, as seen in the comparison between frames (b) and (d). 





7* 

Fig. 3. Semi-quantum Rabi oscillation and photon scattering sequence (upper frames) and corresponding waiting- 
time distribution (lower frames) for four different amplitudes of oscillation in the trap: A =0.0 (a), 0.38 (b), 0.6 
(c), and 0.88 (d). 



A useful way to quantify the behavior as a function of A is through the distribution of times waited 
between two successive scattering events (quantum jumps) — the so-called waiting-time distribution. This 
distribution is plotted for semi-quantum trajectories in the lower series of pictures in Fig. 3, where it is 
narrow (width of order jt = 1) in frames (a) and (c), significantly broadened in frame (b), and broader 
still in frame (d). The distribution depends, most generally, on both the amplitude, A, and phase, 

C = a;Tt-arg(a), (23) 

of the oscillation in the trap at the time of the quantum jump that begins the interval waited r. The 
quantity plotted in Fig. 3 has been averaged over (. For given amplitude and phase, the waiting-time 
distribution is computed within quantum trajectory theory as [39] 

W{r;A,C)=l\c^+\t + T)\\ (24) 



where c^~^\t + r) satisfies Eqs. (9) with initial condition c^~^\t) = 0, \t) = 1. Its mean, 

f(AC)= / drTW{r;A,0, (25) 
Jo 

is sufficient to characterize the metastable states. For semi-quantum trajectories, an approximate 

analytical expression for this mean can be derived as follows. 

The task is to solve Eqs. (9) with the initial internal state the electronic ground state. We first make a 
transformation to consider the complex variable 



+ r) + «c(+) (t + r) = B(r; A, C)e'f^i-'^'0 , 



(26) 



where B{t; A, () and /3(t; A, Q are real functions of r. The equations of motion for Schrodinger 
amplitudes, Eqs. (9), then yield a solution in the form 



B{t; a, C) = exp 

with the phase /3(t; A, Q satisfying the equation 
dp{T;A,0 



[ dr'sin2/3(r';A,C) 
Jo 



(0/2) cos [27rA cos(a;TT + C)] - ^ sin/3(T; A, Q. 



dr 



From Eqs. (24), (26), and (27), the waiting- time distribution is 

W{t;A,C) = -^B^{t;A,0, 
dr 

with mean waiting time, integrating Eq. (25) by parts, 

/•OO 

f(AC)= / dTB'{T;A,0. 
Jo 



(27) 



(28) 



(29) 



(30) 



The development up to this point is exact. Wc now introduce an approximation which holds good 
whenever the mean waiting time is much larger than 7~^. Assuming /3(t; A, Q is very small, such that 
sin/3(T; A, C) « P{t;A,Q, Eq. (27) becomes 



B{t; a, C) = exp 
while Eq. (28) can be solved for 

I3{t;A,0 = (r!/2)Re 



-ll^dT'0\r';A,O 



(31) 



±^ j/2 + inujT 



(32) 



where we have used the standard Fourier series expansion of cos [27r^ cos(wtt + C)] [40] . We then replace 
/3^(t; A, C) in Eq. (31) by its d.c. component (which is independent of the phase C), and from Eqs. (29)- 
(31), the waiting-time distribution is approximated by the exponential distribution 



W{T;A)=f-'e 



_ =-l^-T/f(A) 



with mean waiting time 



7t{A) = 



(0/2)2 f;^2„(2^A) + j2^„(27r^) 



H -1 



n=0 



(7/2)2 + (2„^^)2 



(33) 



(34) 



Note how the dominant term in the series expansion of 7f vanishes for the identified metastable 
amplitudes — when Jq{2ttA) — 0. The approximation (34) is compared with the exact result computed 
from a semi-quantum trajectory average in Fig. 4. The two results are in remarkably good agreement 
everywhere except for amplitudes A < 1 and in a narrow range of amplitudes in between the peaks. 
We see that with increasing center-of-mass energy, the rate of photon scattering on the metastable 
amplitudes very quickly drops by more than two orders of magnitude. Clearly this interplay of the 
coherent oscillation of the ion in the trap and the induced Rabi oscillation must be the principal 
determinant of the heating curves — for small Lamb-Dicke parameters at least — and of the noted 
differences between semi-quantum and quantum models (Fig. 1). 
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Fig. 4. Exact (thick curve) and approximate [Eq. (34)] (tliin curve) semi- 
quantum mean waiting time plotted as a function of the amplitude of 
oscillation in the trap. 

4-2 The semi- quantum diffusion limit: 

Point (i) in Sec. 2 is concerned with the approach of the heating curves for small 77 to a diffusion limit. It 
was noted that the 77-dependence in this limit amounts to a scaling of the time by 77**. We are now in a 
position to elaborate upon these observations. We treat the diffusion limit for semi-quantum trajectories 
first, then compare the description for quantum trajectories in the following section. 

It is convenient to use a set of coordinates rotated to remove the phase of the complex amplitude a. In 
this frame of reference, the momentum "kick" defined by Eq. (lib) is written as defined by Eq. (23)] 



A(0, </), C) = A|| (0, C) + * (0, 0, C) = sin cos 0e'^ , 



(35) 



where A||(0, 0, C) is a "kick" in the amplitude of a and A^ (0,0,0 is a "kick" in its phase. Successive 
"kicks" are distributed randomly, though not necessarily uniformly, over the angles 0, 0, and Q. Taking an 
average over these distributions yields a covariance matrix that depends only on the amplitude of 
oscillation in the trap A: 



C{A) = B(A)B^(A) 



A^(^,0,C) 



A||(0,0,C)A^ (0,0,0 



A||(0,0,OAi (0,0,0 



Ai (0,0,0 



(36) 



The covariance matrix defines the r.m.s. size of the phase-space "kicks" . It is plotted from a simulated 
trajectory average in Fig. 5. To this information we add their average rate — i.e., the inverse of the mean 
waiting time plotted in Fig. 4, also a function only of A. Assuming, then, that the phase-space "kicks", of 
order rj [Eq. (35)], are small, so that significant changes in a accumulate over very many photon scattering 
events, we may put these pieces of information together to arrive at the stochastic differential equation 
with amplitude-dependent diffusion 



1 



B(A) 



dW\\ 



(37) 



where the rotation matrix R[— arg(a)] brings the phase-space increments back to the original coordinate 
system; dW\\ and dW± are independent Wiener increments, with covariances 



dW\idW\i = dW±dW± = dt, dWiidW± = 0. (38) 




The noted 77^-scaling of the time axis in the diffusion hmit [point (i) of Sec. 2] foUows directly from 
Eq. (35)-(37). Note first that B(v4) ^ rj, while, in addition, the oscillation amplitude in wavelengths is 
related to the phase-space amplitude by A — ri\a\/Tr [Eq. (20)]; thus, after multiplying Eq. (37) by rj/n, 
the right-hand side is of order rj^ and the fluctuation variance, at fixed A, grows as rj^t. 

We find that the diffusion model agrees with semi-quantum jump simulations for Lamb-Dicke 
parameters on the order of 77 = 0.2 or less; in fact, there is very little difference for 77 — 0.4 [results for 
rj = 0.2 and 0.4 are indistinguishable in Fig. 1(a)]. Further increase, however, brings noticeably faster 
semi-quantum heating, as shown by the curves for 77 = 0.8, 1.4, 2.2, and 2.6 in Fig. 1(a) [point (iii) of 
Sec. 2]. In this regime the scale of the individual quantum events (momentum "kicks") approaches that 
of the underlying phase-space structure represented by the peaks in Fig. 4. A handful of jumps can 
destabilize the metastable amplitudes much faster than the diffusion process can. 

4-3 Quantum trajectories compared: 

Semi-quantum trajectories are conceptually simple and provide a convenient starting point for our 
understanding. They are nevertheless a rather gross approximation. Points (ii) and (iv) of Sec. 2 record 
two instances where the approximation has a significant effect upon the heating rate. Considering first 
the quantitative difference in the diffusion limit [point (ii)], we might ask how the mean waiting time and 
covariance matrix change as functions of the amplitude of oscillation A. The mean waiting time computed 
from full quantum trajectories is plotted in Fig. 6. For a Lamb-Dicke parameter of 77 = 0.1 it is almost 
indistinguishable from the semi-quantum result. At larger values of 77, however, differences set in: there is 
a gradual smoothing out of the peaks, which is noticeable for 77 = 0.2 and virtually complete for 77 = 1.4. 
This development is explained in Sec. 5. 

Changes to the covariance matrix are more relevant for the diffusion (small 77) limit. Here the most 
important point to note is the difference between Eqs. (11a) and (19a). The former specifies the semi- 
quantum phase-space "kick" , while the latter adds to this the local-frame displacement of the center-of- 
mass wavepacket, (a)iocai, which occurs during the elapsed waiting time following the last quantum jump. 
We might refer to this as the stochastic dipole-force displacement. It arises from the direct action of the 
standing-wave light field on the center-of-mass of the ion. Since resonant excitation constitutes a 
nonperturbative interaction, it is not possible to extract an optical potential to account for this. 
Nevertheless, to add to the momentum "kicks" from fluorescence, there is also momentum exchanged 
between the standing- wave laser field and the ion; this the semi-quantum model omits, while quantum 
trajectories include it as a time evolution of the center-of-mass kets — IV'rec(^)) ^'^'^ I''/'rec(^)) — taking 
place in between the quantum jumps [Eqs. (17)]. 
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Fig. 6. Mean waiting time from quantum trajectories plotted as a function 
of the amplitude of oscillation in the trap; for Lamb-Dickc parameters 
77 = 0.2, 0.4, 0.8 (thinner line to thicker line), and 1.4 (dashed hne). The 
tallest curve (thinest line) is the exact semi-quantum result plotted in 
Fig. 4. 

From Eq. (19b), when characterizing the individual quantum events for quantum trajectories, we 
replace Eq. (35) by 

A(0,0,C,C') = |(a)iocai|e"^' +^??sin^cos</)e<, (39) 

where a dependence enters upon the additional relative phase 

C = arg((a)iocai) - arg(a). (40) 

The covariance matrix is defined as before, with real and imaginary components of A(0, 0, ^, Q') resolved 
in amplitude and phase directions, but now an average over the four phases 9, (j), rj, r]' is taken. The 
numerically computed result is displayed in Fig. 7. Its principal difference when compared to Fig. 5 is that 
the phase variance peaks on the metastable amplitudes to reach a value an order of magnitude larger than 
in the semi-quantum case. Other less obvious changes are observed; for example, the amplitude variance 
also increases to produce a broad background near each metastable amplitude — the likely cause of the 
noted larger heating rate [point (ii) of Sec. 2]. 

The peaking of the phase variance on the metastable amplitudes appears most dramatically when 
phase-space plots are made. Sample plots are shown in Fig. 8, where we plot the complex amplitude a 
from semi-quantum simulations (second column in the figure) and (a)REC from quantum simulations 
(third column). In the diffusion limit (?/ = 0.2), the phase diffusion on the metastable amplitudes 
produces a dramatic "bulls-eye" pattern for quantum trajectories (top right frame), in sharp contrast 
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Fig. 7. Full quantum covariance matrix elements plotted as a function of the amplitude of oscillation in 
the trap: (a) amplitude (thick curve) and phase (thin curve) quadrature variances; (b) amplitude-phase 
cross-correlation. Results computed as a Monte-Carlo average for r; = 0.2; note that a small sampling 
error remains. 
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Fig. 8. Sample semi-quantum (second column) and (niantum (third column) trajectory simulations compared 
through plots of the complex oscillator amplitude in phase space; for Lamb-Dicke parameters 7} = 0.2 (a), 0.8 (b), 
1.4 (c), and 2.2 (d). The first column reproduces the plots from Fig. 2. 



to the semi-quantum result. The phase diffusion is stiU in evidence for r] = 0.8 (second row), but gradually 
ceases to be a feature as r/ is further increased and the diffusion picture breaks down. 

Our final comment on the comparison between semi-quantum and quantum trajectories for small 
Lamb-Dicke parameters is illustrated by Fig. 9. Quantum trajectories provide a quantum description of 
the center-of-mass state — they evolve a center-of-mass wavepacket. What, then, is the form of this state? 
Perhaps we can expect a good approximation to a minimum uncertainty state, which is the closest 
quantum mechanical representation of the phase-space point of Eq. (7). As it turns out, this is not the 
case, even for the smallest Lamb-Dicke parameters considered (though for sufficiently small r] we might 
still expect that result). The states shown as Q functions in Fig. 9 were recovered for a Lamb-Dicke 
parameter rj = 0.2. They illustrate a generally observed amplitude squeezing, with large squeezing in 
the vicinity of the metastable amplitudes and significantly less squeezing at amplitudes in between. The 
degree of squeezing — of the conditional (instantaneously sampled) wavepacket — is plotted as a function of 
the amplitude of oscillation in the trap in the frame to the right in the figure. 

We propose the following quantum measurement interpretation of these results. The ion fluorescence 
may be viewed as the record of an imperfect measurement, one revealing information about the amplitude 
of oscillation in the trap. Most simply, if the fluorescence rate is low, the ion likely oscillates with one 
of the metastable amplitudes; if high, with some amplitude midway in between. Beyond this simplest 
statement, a more quantitative deduction about the actual amplitude of oscillation may be made on the 
basis of the correlation between fluorescence rate and oscillation amplitude given in Figs. 4 and 6. The 
important point to note is that the amplitude discrimination achieved is especially high in the vicinity of 
a metastable amplitude and significantly poorer in between. Where the discrimination is high, the center- 
of-mass wavepacket becomes squeezed so that its predicted uncertainty in the amplitude of oscillation is 
consistent with the data made available though the measurement record (fluorescence)-i.e. the squeezing 




Fig. 9. Squeezing of the centcr-of-mass wavepaclcet as a function of tlic amplitude of oscillation in the trap. The 
scattered photons may be detected and read as the signal of an imperfect measurement of the amplitude A. Near 
a metastable amplitude, the discrimination of the measurement is very good and the centcr-of-mass wavepacket 
becomes squeezed so as to be consistent with the record read (Q function on the lower left): in between the 
metastable amplitudes the discrimination, consequently the squeezing, is not so good (Q function on the upper 
left). To the right, amplitude (thick curve) and phase (thin curve) quadrature variances arc plotted as a function 
of the amplitude of oscillation in the trap. All results are for a Lamb-Dicke parameter r/ = 0.2. 



may be view as the product of a dynamical "collapse of the wavepacket" , which in the case of a near-ideal 
measurement would produce an amplitude eigenstate. 

5. Quantum inhibition of fluorescence and the limit of large quantum jumps 

Points (iii) and (iv) of Sec. 2 note that large Lamb-Dicke parameters yield an ry-dependence in the heating 
curves of Fig. 1 over and above the r/^-scaling of the time axis — the dependence arising from the diffusion 
limit. Point (iv) notes, in particular, that the quantum heating rate decreases for large rj, while the 
inereasc seen for semi-quanttim trajectories seems more reasonable, when, due to their increased size, just 
a few momentum "kicks" are sufficient to "skip over" (destabilize) the metastable amplitudes. In fact, the 
lowered heating rate is caused by a suppression of fluorescence of an entirely different (to the mechanism 
of Sec. 4) quantum mechanical origin, as we now explain. 

Consider, for example, the matrix elements for excitation of the ion out of the ground electronic state 
and the ground state of the trap while absorbing one photon and an even number of "phonons" ; the fi- 
nal state is the excited electronic state and energy eigenstate state of the trap |2n), n = 0, 1, 2, . . .. The 
matrix elements, as a function of t], are 

(2n| cos[kx{t)]\0) = (2n| cos [rj (ae"*'^^* + a^e''"-^*)] |0) 

= (_!)"_/ e"" /^e'^"'^-^*. (41) 
V2n! 

Note that matrix elements for the absorption of an odd number of "phonons" vanish due to the even 
parity of cos(fca;). The squares of matrix elements (41) are plotted for the first few n in Fig. 10(a). From 
the figure we see that transitions |0) |0) dominate in the diffusion limit. On the other hand, for r/ ~ 1, 
the ground-state wavepacket extends beyond ±A/4 where the cosine function changes sign. With the sign 
change the value of the matrix element is reduced. Thus, at 77 = 1.5, the |0) |2) transition is strongest, 
while for ~ 2, |0) ^ |4) and |0) |6) are the strongest transitions. We note then that these strongest. 
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Fig. 10. (a) Matrix elements for transitions out of the ground state of the trap, (0), to the excited state 
2n), plotted as a function of Lamb-Dickc parameter rj for n = 0, 1, 2, and 3 (thicker curve to thinner 
curve); (b) Mean waiting time until the first quantum jump as a function rj. 



multi-phonon transitions are not excited on resonance; the |0) — > \2n) transition is detuned by an 
amount 2nujT — 2nj [Eq. (5)]. The combined effect of these two observations leads to a strong quantum 
suppression of the fluorescence; hence the reduced heating rate of Fig. 1(b). As an iUustration of the 
suppression, in Fig. 10(b) we plot, as a function of ry, the mean waiting time — beginning in the ground 
state of the trap — for the very first photon to be scattered. The waiting time grows rapidly for Lamb- 
Dicke parameters larger than rj = 1.5. 

There is a great deal more to be said about the quantum-stochastic motion of a trapped ion in the 
regime of large Lamb-Dicke parameters. This regime is similar to that of optical frequency cavity QED, 
where the influence of single scattering events is no longer small and the diffusion picture fails, or becomes 
forced at best [32]. As an indication of what is on offer, we conclude with Fig. 11. Here we present 





Fig. 11. Sample states from quantum trajectories for Lamb-Dicke parameters ?7=0.8 (left plots) 1.4 (center plots), 
and 2.2 (right plots). Q functions for the conditional center-of-mass state are plotted. 



Q functions for sample states reached with Lainb-Dicke parameters 7] = 0.8, 1.4, and 2.2. Of course, a 
huge variety of such states are visited throughout the course of a quantum trajectory; the question of what 
interesting quantities might be easily measured has yet to be explored. Note, however, the progression 
from the amphtude squeezed state of Fig. 9 to the states with increasing phase uncertainty of Fig. 11. 
Note also the interesting development of the amplitude uncertainty. The state for ry = 0.8 (on the left in 
the figure), for example, resides on three metastable amplitudes at once, as a superposition. Then there is 
the development from left to right with increasing rj — where, for rj — 2.2, a somewhat disordered phase- 
space structure emerges from the quantum interference pattern. Clearly, there remains much interesting 
behavior to be investigated in future work. 
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